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ABSTRACT 


Extreme  value  distribution  laws  are  obtained  for 
the  lifetimes  of  multi-component  systems  with 
replaceable  components,  under  various  assumptions 
on  the  asymptotic  relationship  between  number  of 
components  in  the  system  and  number  of  spare 
components.  Results  are  given  for  limiting 
distribution  laws  of  order  statistics  from  non- 
homogeneous  samples  and  samples  of  random  size, 
and  applied  to  the  superposition  of  renewal 
processes.  An  attempt  is  made  to  put  extreme 
value  theory  into  a  general  framework  using  the 
notion  of  a  coherent  structure,  and  some  new 
results  utilizing  this  idea  are  presented. 


CONTRIBUTIONS  TO  THE  THEORY  OF  EXTREME  VALUES 


by 

Robert  Harris 

INTRODUCTION 

The  limiting  distribution  of  the  maximum  term  in  a  sequence  of  independ¬ 
ent,  identically  distributed  random  variables  was  completely  analysed  in  a 
series  of  works  by  many  writers,  culminating  in  the  comprehensive  work  of 
Gnedenko  [6].  Results  for  order  statistics  of  fixed  and  increasing  rank  were 
obtained  by  Smirnov  [14],  who  completely  characterized  the  limiting  types  and 
their  domains  of  attraction.  Generalizations  of  these  results  for  the  maximum 
term  have  been  made  by  several  writers;  Juncosa  [10]  dropped  the  assumption  of 
a  common  distribution,  Watson  [15]  proved  that  under  slight  restrictions  the 
limiting  distribution  of  the  maximum  term  in  a  stationary  sequence  of  in¬ 
dependent  random  variables  is  the  same  as  in  the  independent  case,  and  Berman 
[2]  studied  exchangeable  random  variables  and  samples  of  random  size.  A 
bibliography  and  discussion  of  applications  is  contained  in  the  book  by  Gumbel 

[9]. 

This  paper  extends  the  classical  theory.  The  second  section  introduces 
a  model  from  reliability  theory  -  essentially  a  series  system  with  replaceable 
components.  It  is  shown  that  the  asymptotic  distribution  of  system  lifetime 
can  belong  to  one  of  two  types  when  the  number  of  spares  is  fixed  or  of  a 
smaller  order  than  the  total  number  n  of  components,  as  n  becomes  infinite, 
and  that  these  limiting  distributions  are  the  same  as  those  obtained  by 
Gnedenko,  Chibisov  [4]  and  Smirnov.  Sections  3  and  4  deal  wi  :h  the  limiting 
distribution  of  order  statistics  when  the  assumptions  of  common  distribution 
and  fixed  sample  size  are  dropped.  The  results  of  these  three  sections  are 


f 
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then  applied  in  Section  5  to  the  superposition  of  a  large  number  of  renewal 
processes  and  compared  to  the  necessary  and  sufficient  conditions  of 
Grigelionls  [8]  concerning  closeness  to  a  Poisson  process.  In  particular, 
it  is  shown  that  randomness  of  the  sample  size  leads  to  a  point  process  which 
is  a  mixture  of  nonhomogeneous  Poisson  processes.  Finally,  in  Section  6,  an 
attempt  is  made  to  put  extreme  value  theory  into  a  more  general  framework 
using  the  notion  of  a  coherent  structure  [1].  As  noted  there,  most  of  the 
classical  problems  in  extreme  value  theory  are  contained  in  the  following 
general  question:  "Given  a  coherent  structure  with  a  finite  number  of  compon¬ 
ents  and  some  procedure  to  increase  the  number  of  components  without  bound, 
what  are  the  possible  limiting  distributions  for  the  structure  lifetime  for 
given  component  lifetime  distributions?"  For  example,  the  minimum  of  a  set 
of  n  random  variables  corresponds  to  the  lifetime  of  a  series  structure  of 
n  components,  and  the  limiting  procedure  adds  one  component  at  a  time  to  the 
structure.  In  addition,  the  class  of  limiting  distributions  is  characterized 
for  the  case  of  an  arbitrary  coherent  structure  when  the  procedure  for  expan¬ 
sion  is  that  of  repeated  composition. 
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1.  NOTATION  AND  CLASSICAL  RESULTS 

Throughout  this  paper,  the  distribution  function  of  a  random  variable  X 
will  be  denoted  by  P{X  £  x}  ■  F(x)  ,  and  the  tail  of  the  distribution  by. 

P{X  >  x}  -  F(x)  .  The  abbreviation  "d.f."  will  be  used  for  distribution 
function.  A  d.f.  will  be  called  proper  if: 

lim  F(x)  “  1  ,  lim  Fvix)  «=  0 

X-x» 

and  not  all  its  mass  is  concentrated  at  one  point.  Two  d.f.'s  F^(x)  a.d 

?2(x)  are  said  to  be  of  the  same  type  if  there  exist  constants  A  >  0  and  B 

such  that:  F^(Ax  +  B)  ■  F£(x)  for  all  values  of  x  .  Unless  otherwise  stated, 

all  d.f.'s  will  be  assumed  proper  and  all  limiting  d.f.'s  should  be  taken  to 

mean  limiting  types  of  d.f.'s.  Let  X^,X2»  ....  X^,  ...  be  a  sequence  of 

Independent  random  variables  with  common  distribution  F(x)  ,  and  let 

Cn  “  minCX^.Xj . X^)  .  Then  the  limiting  d.f.  of  belongs  to  exactly 

one  of  three  types,  [6];  that  is  to  say,  if  there  exist  sequences  of  normalizing 

constants  {a  >  0}  and  {b  }  and  a  d.f.  G(x)  such  that: 
n  n 

lim  pja  1  (£  -  b  )  <  x  |  ■  G (x) 

n-*»  '  ' 

at  each  continuity  point  of  G(x)  ,  then  G(x)  belongs  to  one  of  the  following 
types: 


\l)M 

-  0 

for 

x 

< 

0 

»  1-exp  C-x°] 

for 

X 

> 

0  ,  a  >  0 

(1.1) 

*(2)(x) 

-  1-exp  [-(-x)'a] 

for 

X 

< 

0  ,  a  >  0 

-  1 

for 

X 

> 

0 

♦  ^(x)  -  1-exp  [-exp  x] 


-HO  <  x  <  00 
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The  domain  of  attraction  ui  a  limiting  d.f.  G(x)  is  the  set  of  all  d.f.'s 

F(x)  such  that  for  suitable  choice  of  normalizing  constants  {a^  >  0} 

and  {b  } 
n 

(1.2)  lim  F°(a  x  +  b  )  ■  G(x)  . 

n  n 

n-H» 

By  a  well-known  theorem  of  Khintchine  (e.g.,  see  [7],  p.  40),  each  d.f.  can 
belong  to  at  most  one  domain  of  attraction.  Necessary  and  sufficient  condi¬ 
tions  were  given  by  Gnedenko  [6]  for  a  d.f.  to  belong  to  the  domain  of 
attraction  of  *(2)^  or  *(3)^  *  For  example>  F(x)  is  i° 

the  domain  of  attraction  of  4>^(x)  if  and  only  if  3  x0  such  that 
F(xq)  ■  0  ,  F(xo  +  e)  >  0  for  each  c  >  0  and 

F(x  +  tx) 

(1.3)  lin  -  t°  for  all  t  >  0  . 

x->0+  '>xo  +  ' 

The  k**1  smallest  variable  from  (X^^,  . ..,  X  )  will  be  denoted  by  % 

so  that  £  ^  ■  £  ;  limiting  d.f.‘s  for  these  random  variables  as  obtained 

n  n 

by  Smirnov  and  Chibisov  will  be  introduced  as  needed. 

The  notation  f(x)  -  0(g(x))  as  x  •*  a  will  mean  |f(x)/g(x)|  k  <  *> 
as  x  a  ,  and  f(x)  ■  o(g (>*))  as  x  ■+•  a  will  be  used  to  denote  that 
f  (x)/g(x)  -*■  0  as  x  -*•  a  .  Likewise  f(x)  -  g(x)  as  x  -►  a  implies  that 

f  (x)/g(x)  1  as  x  a  . 


* 


* 


♦  (x)  will  be  used  for  the  standard  normal  (0,1)  d.f. 


... 
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2.  STRUCTURES  WITH  REPLACEMENT 

The  problem  that  is  investigated  here  is  the  following:  a  system  con¬ 
sists  of  n  identical  and  independent  components  in  series,  with  in  inactive 
spare  components  available  which  instantaneously  replace  the  components  as 
they  fail,  until  there  are  no  more  spares,  whereupon  the  system  falls.  The 
system  lifetime  will  be  denoted  by  ,  (m  +  1)  being  the  total  number 

of  component  failures  which  must  occur  before  system  failure.  The  investiga¬ 
tion  is  in  two  parts,  corresponding  to  the  cases  when  m  -  m(n)  is  of  a 
smaller  order  than  n  or  of  the  same  order  as  n  ,  and  a  third  subsection 
describes  how  some  of  the  results  may  be  carried  over  to  more  general  types 
of  systems.  It  is  assumed  in  this  section  ti  :*.  F(O-)  «  0  . 

Extreme  Terms 

Let  G  (x)  *  pL(m+l)  <  xj  .  Then  it  is  shown  that  the  class  of  limiting 
nm  (  n  —  ) 

d.f.'s  for  the  system  lifetime  as  n  °»  ,  with  appropriate  linear  norming 
constants,  is  the  same  as  the  limiting  d.f.’s  of  the  corresponding  order 
statistics  provided  that  m  is  finite  or  of  smaller  order  than  Sn  ,  as  in 
the  following  two  theorems. 


Theorem  2.1: 

* 

The  limit  laws  for  sequences  Gntn(®nx  +  ^  )  of  system  lifetime  d.f.'s, 
with  m  fixed,  are  exhausted  by  the  following  two  types: 


(2.1) 


(m 


1  f  -v  m-1 

“ Ti  J  e  y 


dy 


*(3)  W 


(m 


w 


-y  m-1  , 
e  7  y  dy 


for  x  <  0 


for  x  >  0  ,  a  >  0 


00  <  X  <  00 
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Theorem  2.2; 

If  b  “  cna  ,  with  c>0,0<a<!s,  then  the  only  possible  limit 

h 

d.f.  s  for  the  sequence  G  (a  x  +  b  )  are: 

n  nm  n  n 


G(d  (x)  -  4>(x) 

(2.2) 

G^2) (x)  "  0 

for  x  0 

-  4>(Blogx) 

for  x  >  0  , 

Notice  that  G^(x) 

is  the  log  normal  d.f. 

Some  preliminary  results  are  needed  before  the  proofs  of  Theorems  2.1 
and  2.2  can  be  given: 


(2.3) 


G*  (x) 
nm 


m 

l 

J-o 


*1  + 


n 

n 

+  in«j  k-1 


(i  )  (i  +1) 

F  *  (x)  -  F  *  (x) 


(k) 

Where  F  (x)  is  the  k-fold  convolution  of  the  d.f.  F(x)  and  the  inner 

summation  is  over  all  nonnegative  combinations  of  (1. ,1.,  ....  i  )  which 

1  l  n 

sum  to  J  .  This  formula  follows  from  the  superposition  of  n  identical 
renewal  processes. 

The  d.f.  F(x)  will  be  assumed  to  be  concentrated  on  the  nonnegative 
real  axis  in  this  section  since  the  concept  of  component  lifetime  is  meaning¬ 
ful  only  in  this  case.  Use  will  be  made  of  the  Inequality 


(2. A)  F(k)(x)  <  (F(x)}k  ,  V  k>l,  V  x^O. 

It  is  convenient  to  speak  of  n  "sockets"  in  series,  each  of  which  must 
contain  a  working  component  for  the  system  to  work.  When  m  is  not  too 
large,  a  key  step  in  the  proofs  will  be  to  show  that  the  probability  of  two 
or  more  failures  in  any  socket  is  negligible  as  n  -►  ®  .  Define 
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(2.5)  G  (x)  ■=  l  (")  Fn"J(x)  (x) 

nm  j=o  Xil 

i.e.  ,  the  survival  probability  of  an  (m  +  1)  -out-of-n  system. 

Theorem  2.3: 

L 

If  ra  -  o(n  )  as  n  -*•  °°  ,  and  if  {a  >  0}  and  {b  }  are  sequences  of 

n  n 

normalizing  constants  such  that 

-k 

F(a  x  +  b  )  *=  o(n  )  as  n  -*  «  ,  V  x  >  0  ,  t.;en 
n  n  *  —  ’ 

(2.6)  lim  Ig  (a  x  +  b  )  -  G  (a  x  +  b  )  I  «  0  ,  V  x>0. 

1  nm  n  n  nm  n  n  1  ’  — 

n-x» 

The  proof  of  this  theorem  will  depend  on  the  following  lemmas. 

Lemma  2.1: 

(i)  The  number  of  ways  in  which  j  failures  can  occur,  in  such  a  way 
that  at  most  one  failure  occurs  in  each  socket,  is  (j)  • 

(ii)  The  total  number  of  ways  in  which  j  failures  can  occur,  the 
number  of  failures  in  any  socket  being  arbitrary,  is 

Proof ; 

(i)  Follows  by  considering  the  coefficient  of  in  (1  +  Z)n  . 

(ii)  Follows  by  considering  the  coefficient  of  Z^  in 
(1  +  Z  +  Z2  +  ...)n  =  (1  -  z)"n  . 

Lemma  2.2: 

If  0  1  j  £  b  ,  and  m  -  o(n  )  as  n  -*■  °°  ,  then 


("T1)- 


( j )/(n  +  j  '  l)  *  1  ■  as  n  *  •  • 
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Proof : 


• £  Gir  n- vr 


1  (‘  -  5)' 

*  (■  -if  ■ 


for  n  sufficiently  large. 


-  exp.  l-S)  , 

where  e  >  0  is  arbitrary.  The  result  follows  on  letting 
Proof  of  Theorem  2.1; 

Deff-ne  the  following  notation: 

1)  A*-(n+ri)-Q 

2)  unj  (x)  -  Fn"j  (x)  (x) 

3)  vnj(x)  -  Fn~J(x)  {F(x)  -  F(2)(x)}J 

4)  Wnj(x)  w111  be  used  for  all  terras  of  the  form: 


e  ■*  0  . 


n  (  (V  (m-d  ) 

^  fr  M  -  F  k  Mj  , 


Where  ij  +  . . .  +  ln  .  <  ,„d  at  least  one  of  the  1  >  2  .  Notloe 


from  (2.4)  that: 


(2.7) 


°-Vx)  iunj(*)  •  »1«,JW<.  (x)  . 


G*  (X)/G„  (x) 
nn  nra 
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1 

! 


But : 


Also : 


m  i  m 

<  l  A  w  (x)  /  j  (")u  .(x 
-  joO  n*  nJ  /j-o  nJ 


i  n  Vn3(X>/(")UnJ<’l) 

J  »  »  »  1,1 


-»■  0  ,  by  (2.7)  and  Lemma  2.2. 


°il_  joO"3(l,)/jo 

■J!0(j)IUn3<’,)'V"J<X)1/j0^Uni<X) 

.<  max  [u  (x)  ~  vni  ^x^^uni 

-  1  -  [1  -  F(2)(x)/F(x))m 

<_  1  -  [1  -  F(x)]m  ,  by  use  of  (2.4)  . 


✓ 


Now  if  x  is  replaced  by  (a^x  +  bn)  and  the  second  assumption  of  the 
theorem  used,  it  is  seen  that  the  last  term  approaches  zero  as  n  ■+  00  . 
Combining  results: 


(a  x  +  b  )/G  (a  x 
n  n  nm  n 


+  v  - 1| 


0  ,  as  n  -*•  *  , 


and  since  d.f.'s  are  bounded,  the  theorem  is  proved. 

Proof  of  Theorems  2.1  and  2.2: 

To  examine  the  possible  normalizing  sequences  {a^}  and  (b^)  which 
satisfy  the  conditions  of  Theorem  2.3,  it  is  necessary  to  consider  separately 
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the  cases  where  in  remains  finite  or  m  -*•  »  .  Suppose  first  that  m  remains 
finite.  Then  Smirnov  [1A]  has  shown  that  in  order  for 

(2.8)  G  (a  x  +  b  )  -*•  G(x) 

nm  n  n 

for  suitable  choice  of  normalizing  constants,  where  G(x)  is  a  proper  d.f., 
it  is  necessary  and  sufficient  that 

(2.9)  v  (x)  -  nF(a  +  b  )  -*■  v(x) 

n  n  n 

where  v(x)  is  a  nondecreasing  nonnegative  function  defined  by: 

v(x) 

(2.10)  (m  -7)7  /  e  y  ym  1  dy  **  G(x) 

*/o 

Furthermore,  he  proved  that  (up  to  a  linear  transformation)  the  function 
v(x)  must  be  one  of  the  three  forms  x°  ,  (-x)  a  or  e”X  ,  where  a  is  an 
arbitrary  positive  constant.  The  domain  of  attraction  corresponding  to  the 
second  form  for  v(x)  consists  of  d.f.'s  which  are  unbounded  below,  so  that 
on  using  Theorem  2.3  and  the  nonnegativity  assumption  on  the  {X^}  ,  Theorem 
2.1  is  proved. 

Now  suppose  that  m  -  cna  ,  where  c>0,0<a<  Chibisov  [A] 

has  shown  that  G  (a  x  +  b  )  -*■  G(x)  if  and  only  if 
nm  n  n  J 

(2.11)  u  (x)  «  m  **  {nF(a  x  +  b  )  -  m)  -*•  u(x) 

n  n  n 

where  u(x)  is  defined  by  the  equation. 

(2.12)  G(x)  -  *(u(x)) 


and  $  is  the  normal  (0,1)  d.f.  The  function  u(x)  must  be  of  the  same 
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type  as  one  of  x  ,  glogx  or  -01og|x|  ,  where  6  >  0  is  an  arbitrary 
constant,  and  the  domain  of  attraction  corresponding  to  the  third  form  contains 
only  d.f.'s  which  are  unbounded  below.  For  a  normalizing  sequence  which  sat¬ 
isfies  (2.11),  it  is  clear  that  FCa^x  +  b^)  ■  0(na  -  o(n  **)  ;  thus  the 
conditions  of  Theorem  2.3  are  satisfied  and  Theorem  2.2  is  proved. 

Similarly,  characterisations  of  the  domains  of  attraction  of  these  limit 
d.f.'s  may  be  made.  Note  also  that  one  might  wish  to  restrict  the  limiting 
1 'w  itself  to  correspond  to  a  nonnegative  random  variable,  thus  eliminating 
one  of  the  types  in  Theorems  2.1  and  2.2. 

The  assumption  that  the  spares  have  the  same  lifetime  d.f.  as  the  original 

A  A  _ .  L 

components  is  unnecessary;  any  d.f.  F  (x)  such  that  F  (anx  +  bn>  *=  o(n 

will  suffice.  The  appropriate  modifications  to  the  proof  of  Theorem  2.3 

present  no  difficulty. 

It  would  be  desirable  to  relax  the  restriction  a  <  %  which  appears  in 
the  conditions  of  Theorem  2.2.  Results  may  be  obtained  for  a  <  2/3  as 
described  by  Lemmas  2.3  and  2. A,  but  the  more  general  case  a  <  1  does  not 
seem  amenable  to  analysis  and  a  counter-intuitive  reason  for  this  is  given  in 
Lemmas  2.5  and  2.6. 

"st" 

Let  the  symbol  _>  stand  for  "stochastically  greater  than". 

Lemma  2,3: 

For  independent,  identically  distributed  normegative  component  random 
variables : 

f(m)  st  Cm)  st  . (m) 

5n  -  "n  -  5n+m 


% 


where  the  £'s  are  the  corresponding  order  statistics. 


> 


I 
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Proof : 

The  first  part  of  the  inequality  follows  by  observing  the  replacements 
themselves  may  fail,  thus  giving  rise  to  more  failures;  the  second  part  by  .  , 

observing  that  time  to  system  failure  decreases  if  the  spares  are  subject  to 
failure  from  the  initial  instant. 

Lemma  2 . A : 

If  m  -  cn°  ,  with  c  >  0  ,  ^  a  <  2/3  ,  then  the  limit  d.f.'s  (2.2) 

* 

are  possible  for  the  sequence  G^Ca^x  +  bn)  •  I 

Proof: 

Suppose  that  F  and  {a  >0},  {b  }  are  such  that  (2.11)  and  (2.12)  hold, 

n  n 

bo  that  I 

i 

i 

•  I 

(2.13)  F(a  x  +  b  )  “=  m/n  +  u(x)/m/n  +  o^/n)  . 

n  n 

Then 

(n  +  m)  F(a  x  +  b  )  -  m 

(2.14)  -  -  u(x)  +  0(m  '  /n)  ■>  u(x)  . 

/£  ; 

Thus  both  an^  -  bn)  and  -  b  )  have  the  same  limiting  d.f. 

and  hence  by  Lemma  2.3  so  does  a  ^  (o^01^  -  b  )  .  Thus  with  Chibisov's  results, 
the  lemma  is  proved. 

It  should  be  noted  that  although  Lemma  2.4  shows  that  the  limiting  d.f.'s 
(2.2)  are  possible,  it  does  not  rule  out  other  limiting  d.f.'s,  in  constrast 
to  the  results  of  Theorems  2.1  and  2.2. 

w* 

Lemma  2.5: 

The  number  of  ways  that  m  failures  .an  occur  in  n  sockets  with  at 

I 

most  r  failures  per  socket  is 


f  ifiWC'*,, 


where  [x]  denotes  the  largest  integer  less  than  or  equal  to  x  . 


Troof : 


The  form  of  c(n,m,r)  follows  by  observing  that  it  is  the  coefficient 


of  Z  in: 


(1  +  Z  +  Z2  +  ...  +  Zr)n  -  (1  -  zr+1)n  (1  -  Z)"n  . 


Let  c(n,m)  = 


(B+:‘1)  -- 


the  total  number  of  ways  that  m  failures 


can  occur  in  n  sockets. 


Lemma  2.6: 


If  m  ~  cn  ,  where  c>0,0<a<l,  and  r  is  fixed,  then 
c(n,m,r)/c(n,m)  -*■  1  as  n  -*■  ®  provided  r  >_  [1/(1  -  a)]  . 


Proof : 


Write  c(n,m,r)  =  aQ  - 
and  s  -  [m/(r  +1)]  .  Then: 


<->e  s'  •  -*•  S  ■  GKT-1) 


1  _  1  +  1  _  (n+m-ri-i-1) !  (m-ri-r-i-2) ! 

...  n  -  i  (m-ri-i-1)!  ’  (n-m-ri-r-i-2) ! 


>  L±A  h  + - s — 

-  n  -  i  m-ri-i-1 


i  t  r+1  , 

>  A  1  +  -S-r  -  -/r  ’  “y  ' 

—  n  (  o-l  p  (n) 


Then  p  (n)  ■+•  0  as  n  -*•  »  provided  (r  +  1)  >  (1  -  a)  .  Wow 
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|c(n,m,r)/c(n,m)  -  lj  «=  \bQ  “  ^  +  a2  “  •••  (-)S«8)/a0  "  1 ' 

<  I (a.  +  a0  +  . . .  +  a  )/a  | 

—  1  l  so 

2  g 

_<_  p  +  p  +  . . .  +  p 

-  p(l  -  ps)/(l  -  p)  ►  0  as  n  -*■  ®  . 

Thus  for  the  case  m  -  cn°  ,  0  <  a  <  1  ,  although  It  is  tempting  to  think 
that  one  need  only  consider  at  most  o.ie  failure  per  socket  in  the  limit, 

Lemma  2.6  shows  that  a  large  number  of  the  total  ways  of  failure  actually  in¬ 
volve  more  than  one  failure  per  socket. 

Central  Terms 

The  results  obtained  in  the  first  part  of  this  section  are  for  the 
limiting  d.f.'s  of  extreme  terms  in  which  the  number  of  spares  is  of  a  smaller 
order  than  the  number  of  components  in  the  system;  this  part  treats  the  cen¬ 
tral  terms  where  the  numbers  of  spares  and  components  are  of  the  same  order. 

It  is  shown  in  Theorem  2.5  that  under  fairly  weak  conditions  the  limiting 
d.f.  of  a  *  (n ^  -  b  )  ,  for  appropriate  choice  of  a  >  0  and  b  ,  is 
the  normal  d.f.  For  simplicity  of  notation,  it  is  assumed  that  m  *»  n  -  1 
although  it  is  obvious  that  Theorems  2.4  and  2.5  hold  with  slight  modifica¬ 
tions  when  m  -  m(n)  is  such  that  m(n)/n  -*■  X  ,  0  <  X  <  «  . 

Definition: 

Following  Kolmogorov  [11]  and  Smirnov  [14],  a  sequence  {X^}  of  random 

variables  is  said  to  be  stable  if  3  constants  a  such  that 

n 

p{|Xn  -  a  |  <  e}  -  1  ,  as  n  -*■  °°  ,  for  each  fixed  e  >  0  . 

Theorem  2.4  demonstrates  the  stability  of  the  sequence  of  system  life¬ 
time  jn^nM  unc*er  mild  restrictions.  Some  additional  notation  is  needed; 


let  M^(t)  denote  the  number  of  component  failures  in  the  socket  up  to 

n 

and  including  time  t  ,  1  _<  i  <  n  ,  S  (t)  -  £  N  (t)  the  total  number  of 

n  i-1  1 

2 

failures.  Set  y(t)  -  E{N^(t)}  and  a  (t)  ■  Var{N^(t))  as  the  mean  and. 
variance  of  N^(t)  .  It  is  veil  known  that  renewal  counting  functions  N^(t) 
have  finite  moments  of  all  orders  for  each  fixed  t  so  that  the  existence 
of  w(t)  and  o(t)  is  guaranteed. 

Theorem  2. A; 

If  p(t)  is  increasing  in  a  neighborhood  of  t  ■  g  *  (1)  ,  then  the 
sequence  | |  is  stable. 

Proof: 

Fix  c  >  0  and  let  t^  be  the  unique  t  such  that  p(t)  -  1  .  Then 
jtj^  -  c  <  n*”*  <  t1  +  c  <->  [sR(t1  +  c)/n  >  1  >  sn(tx  -  c)/nj  . 

For  arbitrary  c  >  0  ,  P{|SQ(t)/n  -  p(t)|  <  e*}  1  ,  as  n  •  ,  for  all 

finite  t  ,  by  the  weak  law  of  large  numbers.  Thus 

P{Sn(ti  +  e)/n  >  u(t^  +  c)  -  e  }  *►  1  ,  and  by  choosing  c  sufficiently  small, 

it  is  clear  that  w(t.  +  e)  -  c  >1  and  so  P{S  (t.  +  e)/n  >  1}  ■*  1  . 

l  ~  n  x 

Similarly,  P{Sn(t^  -  c)/n  <  1}  -*■  1  ,  so  that  finally 

(2.15)  P|lnnn)  '  hi  <  f-|  -  1  • 

In  xact,  Theorem  2. A  can  be  replaced  by  a  stronger  result  that  is 
analogous  to  the  strong  law  of  large  numbers,  viz.  pjlim  ■  tjJ  ■  1  . 

The  proof  of  this  is  similar  to  that  of  Theorem  2. A  with  the  strong  lav  of 
large  numbers  applied  to  the  sum  Sn<t)  . 


•  -*»»»** 
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Theorem  2.5; 

If  y(t)  has  a  positive  first  derivative  y'(t)  at  t^  then 
/n(n^n)  -  tj  has  a  limiting  normal  d.f.  with  mean  zero  and  variance 
Mt^/v'Up}2  . 

The  proof  of  this  theorem  depends  on: 


Lemma  2.7: 


If  y(t)  is  continuous  at  some  point  t^  ,  then  o(t)  is  continuous 


at  t^  . 


Proof: 

Assume  that  y(t)  is  continuous  at  t^  ;  then  from  the  following 
representation  (e.g.,  see  [1],  p.  54) 


(2.16) 


o2(t) 


f 


2  I  y(t 
'0 


-  x)dy(x)  +  y(t)  -  (y (t) }“ 


it  suffices  to  show  that  the  first  term  in  this  expression  —  denoted  by 
v(t)  —  is  continuous  at  .  For  small  h 


(2.17) 


v(tj^  +  h)  -  v(tj^)  » 


r  } 

m  I  w(tj+h-x)dy(x)  +  /  [ p  (t j+h-x)  -  y(t^  - 


x) ]dy (x)  . 


The  first  term  in  (2.16)  is  not  greater  than  y(h)[y(t^  +  h)  -  y(t^)]  , 
which  becomes  arbitrarily  small  as  h  0  by  continuity  of  y  at  the  point 
.  The  second  term  will  give  a  nonnegligible  contribution  only  if  y  has 
discontinuities  at  both  t^  -  x  and  x  ,  for  some  x  ,  0  <  x  <  t^  .  However, 


since 
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(2.18)  y(t)  -  l  F(1)(t)  , 

1-1 

this  means  that  F^(t)  is  discontinuous  at  t^  -  x  ,  for  some  i  ,  and. 
F^*0(t)  discontinuous  at  x  ,  for  some  j  ,  so  that  F^+^(t)  is  dis¬ 
continuous  at  t^  and  this  with  (2.18)  contradicts  the  hypothesis.  Thus, 
a  is  continuous  at  the  point  t^  . 

Proof  of  Theorem  2.5: 

For  fixed  x 


(2.19) 

o> 


S  (t.  +  x/ /n)  -  ny(t  +  x//n) 
n  x  l _ 

/no(tjL  +  x//n) 


n  -  ny(t^  +  x//n) 

^no(t.  +  x//n) 

1 


n 

Now  S  (t1  +  x//n)  may  be  written  in  the  form  J  X  ,  ,  where 

n  i.  •  ■  nK 

k*=l 

e  N^(t^  +  x/i/n)  ;  it  is  clear  that  the  {X^}  are  independent, 
identically  distributed  and  have  finite  moments  of  all  orders.  Thus,  a 
modification  of  Liapunov's  version  of  the  central  limit  theorem 
(see  [12],  p.  277)  may  be  applied  to  give: 


(2.20) 


^(^  +  x/ /n)  -  ny^  +  x/v^ri) 
i/na(t^  +  x/i^n) 


<  u 


4>(u) 


where,  as  before,  $  is  the  normal  (0,1)  d.f. 

Now  y(t^  +  x//n)  may  be  written  in  the  form 

y(t1  +  x//n)  -  y(t1)  +  (x//n)y' (t^  +  oQ/Zn)  , 


(2.21) 


as  n  -*  00  .  Also,  from  Lemma  2.7, 


(2.22) 


' (t  1  +  x/t/n)  -*•  o^)  . 


Combining  (2.21)  and  (2.22) 


(2.23) 


n  -  nu(t1  +  x/v^n)  /xpMtjA 
/no ( t ,  +  x/v^n)  \  °^ti^  / 


Since  the  normal  d.f.  is  continuous,  the  conclusion  of  the  theorem  follows  by 
substituting  (2.20)  and  (2.23)  into  (2.19). 

It  should  be  noted  that  the  proof  of  the  theorem  is  not  sensitive  to  the 
assumption  of  common  lifetime  d.f.  for  r^ch  of  the  original  and  spare  components. 
All  that  is  needed  is  a  central  limit  theorem  to  hold  for  the  sum 
Sn(t^  +  x//n)  and  convergev.ce  of  the  appropriate  sequence  of  constants  as 
in  (2.23). 


Examples: 

-Xt  2 

1.  Suppose  that  F(t)  *  1  -  e  ,  so  that  y (t)  *■  a  (t)  *  Xt  .  Then 


/"(nn 


(n)  .  -IV 
n  > 


has  a  limiting  normal  d.f.  with  mean  zero  and 


variance  X  .  In  fact,  this  result  can  be  obtained  quite  simply 
by  observing  that  the  times  between  consecutive  failures  are  in¬ 
dependent,  identically  distributed  exponential  random  variables. 

2.  Nonider.tical  components.  Suppose  that  the  original  components  have 


lifetime  d.f.  Fg(t)  ■  Xy  F(x)dx  and  the  spares  have  lifetime  d.f. 

-2Xt 

F(t)  ,  where  F(t)  ■>  1  -  (1  +  2Xt)e  is  a  gamma  d.f.  Thus,  the 

sequence  of  failures  in  each  socket  corresponds  to  an  equilibrium 
renewal  process,  so  that  w(t)  *=  Xt  and  o^(t)  ■>  Xt/2  +  1/8  . 
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Then  /n  (n^n)  -  A-1) 
variance  (1/2  +•  1/8  e 


has  a  limiting  normal  d.f.  with  mean  zero  and 


k-out-of-n  Structures 

The  methods  of  this  section  c..u  be  applied  to  more  general  types  of 
systems  with  replaceable  components.  For  example,  consider  a  k-out-of-n 
system  with  m  spares  where,  as  before,  the  component  lifetimes  are  assumed 
to  be  independent  and  identically  distributed.  As  components  fail,  they  are 
immediately  detected  and  replaced  by  new  components  until  m  replacements 
have  been  made;  the  system  fails  when  k  additional  failures  have  occurred, 
i.e.,  k  +  m  in  all.  Let  the  system  lifetime  be  denoted  by  . 

Lemma  2.8: 

If  k  +  m  <  n  ,  then 

(k+m)  $t  (ra)  <^t  (k) 

^n  —  ^nk  —  St 

Proof : 

The  first  part  of  the  inequality  is  proved  as  in  Lemma  2.3,  and  the 
second  part  by  observing  that  a  system  with  spares  survives  longer  than  a 
system  without  spares. 


Lemma  2.9: 

If  k  +-  m  <  n  ,  then 


. (k+m)  st  „  (m)  st  (k+m) 
’n  -  Cnk  -  nn 


% 
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Proof : 

Only  the  second  part  remains  to  be  proved;  this  follows  by  noting  that 

between  the  m  and  (k  +  m)  failures  there  are  fewe:  than  n  components 

liable  to  failure  and  so  system  failure  is  stochastically  larger  than  in  the 
case  where  replacements  are  continually  available. 

Making  certain  assumptions  about  the  behaviour  of  k  -  k(n)  and  m  =  m(n) 
as  n  -*•  *  enables  some  deductions  to  be  made  concerning  the  limiting  d.f.'s 
of  an^  -  bn )  .  For  example,  consider  the  two  cases: 

(i)  m/n  -*  0  ,  k/n  ->  A  ,  0  <  A  <  1  .  By  using  Lemma  2.8  and  the  results 
of  Smirnov  [ 14 ]  concerning  limit  d.f.’s  of  central  order  statistics, 
the  limiting  d.f.'s  of  system  lifetime  may  be  completely  character¬ 
ised. 

(ii)  (k  +  m)  finite  or  (k  +  m)  ~  cna  ,  c  >  0  ,  0  <  a  <  1/2  .  Then 

Lemma  2.9  and  Theorems  2.1  and  2.2  enable  one  to  describe  completely 

the  possible  limiting  d.f.’s. 
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3.  ORDER  STATISTICS  FROM  A  NONHOMOCENEOUS  SAMPLE 

In  this  section,  it  is  shown  that  the  limiting  d.f.'s  for  order  statistics 
of  fixed  rank  from  an  independent  but  not  identically  distributed  sample  are 
essentially  derived  from  the  limiting  d.f.  of  the  minimum  as  in  the  identical 
case.  Suppose  that  has  d.f.  F^(x)  »  then  Juncosa  [10]  has  shown 


Lemma  3.1: 


If  A(x)  is  a  positive,  nonconstant,  nondecreasing  function  such  that 

for  suitable  {a  >  0}  and  (b  }  : 

n  n 


n 


(3.1) 

and 


lim  £  F,  (a  x  +  b  )  «  A(x) 

km  k  n  n 
“1 


max  F.  (a  x  +  b  )  +  0  ,  V  x, 
,  .  k  n  n  ’  ’ 

l<k<n 


then 


n 

lim  P{£  >  a  x  +  b  }  **  lim  II  F,  (a  x  +  b  )  ■  exp[-A(x)]  . 

nn  n  ,  ,  k  n  n 

k=l 


Theorem  3.1: 

If  the  conditions  of  Lemma  3.1  are  satisfied,  and  m  is  fixed,  then 


The  following  lemmas  are  needed : 

Lemma  3.2: 

>  0  ,  v,  >0  are  such  that: 
—  kn  — 


lim  P 
n-*» 


.  (m) 


a  x 
n 


+  b 


(x) 

jl 


|  exp [-A (x)  ]  . 


% 


If 


t 
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(1)  Ukn  "  U  *  uniformly  in  k  ,  l<k<n  ,  and 
n 

(ii)  ^  vkn  ^  v  *  as  n  "■  “  »  where  u  and  v  are  finite,  then: 

k®l 


(3.2) 


J.Uknvkn  "  uv  •  as  n 
k=l 


Lemma  3. 3: 


If 


d  >0 
kn  — 


n  ,  then: 


is  such  that  max  d 
l<k<n 


kn 


and  £  d 


k=l 


kn 


A  (x) 


as 


(3.3) 


lim  l  \  „  •  •  •  d,.  _  -  4^ 


k^n  k^n 


y 


i 

where  £  is  a  sum  over  all  k±  such  that  1  1  ^  <  k£  <  .  .  .  <  k  <  n  . 


Proof : 


The  result  holds  for  J  -  1  ,  so  assume  the  inductive  hypothesis  that  it 
holds  up  to  J  -  1  .  Clearly 


(3.4) 


LWJ  *  aJ  00  ■ 


Also 


u 


k-1 


kn  j 


(  Terms  containing 


y  I  euuLaining 

3  1  VnVn  dk  n  +  4  j  squares  and  higher 

1  1  J  /  powers  of  d, 

V  kn 


It  is  easy  to  see  that  the  second  term  above  is  asymptotically  negligible. 
For  example: 


^  dkindk2n  * 


d.  <  max  d,  •  V  d  d  A 

kj-i""i<k<„  kn  1  y  y  kj-i" 
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and  Che  right  hand  side  tends  to  zero  from  the  inductive  hypothesis  and  the 
assumptions  of  the  lemma.  The  conclusion  follows  by  comparison  of  3.4  and 


Proof  of  Theorem  3.1: 


(3.6) 


,<x)  "  P!i;nm)  ’  *!  "  J0Adn<x)  • 


where 


(3.7) 


"  A, 


A  (x)  -  l  n  F  (x)  F  (x)  , 


v* 

and  )  is  the  summation  over  all  terms  such  that  i,  +  i_  +  ...  +  i  *  i 
L  1  2  n  J 

r* 

and  i.  »  0,  1  .  Now  A,  (a  x  +  b  )  may  be  written  in  the  form  )  u  v 

k  Jn  n  n  J  L  rn  rn 


where 


(3.8) 


n  1-i,  n  i, 

a  -  n  F.  ,  v  -  n  F,  , 

rn  .  .  k  rn  .  ,  k 

k=l  k=l 


r  is  a  symbol  for  the  partition  r  «=  (i^i^ . i  )  >  and  the  argument  is 

emitted  but  understood  to  be  a  x  +  b  .  Since  u  differs  from 

n  n  rn 

n 

II  x  +  b  )  by  only  a  finite  number  of  factors,  Lemma  3.1  implies 

k-1 


(3.9) 


u  (a  x  +  b  )  -*•  exp.  [-A(x)]  , 
rn  n  n  r  ’ 


uniformly  in  r  . 


If  d^n  is  identified  with  F^Ca^x  +  b  )  ,  the  conditions  of  Lemma  3.3  are 


satisfied,  and  so 


(3.10) 


V*  /  ,  ,  s  (x) 

)  v  (a  x  +  b  )  -*■  —rr  • 
L  rn  n  n  j 


2  k 


Application  of  Lemma  3.2  together  with  3.9  and  3.10  then  gives 
(3.H)  Ajn(V  +  V  _>  exPt"A(x>)  •  • 


The  proof  of  the  theorem  is  complete. 
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4.  ORDER  STATISTICS  FROM  A  SAMPLE  OF  RANDOM  SIZE 

Consider  a  sample  ‘ ‘ '  XN  j  independent  random  variables 

with  common  d.f.  F(x)  ,  where  N  is  a  random  variable  distributed 

n 

independently  of  the  {X^}  .  Assuming  that  Nn  ■*  03  in  probability  as 
n  •*  ®  ,  the  order  statistics  are  defined  with  probability  one  for 

fixed  m  as  n  -*■  00  .  To  avoid  writing  conditional  expectations  it  is  assumed 

below,  without  loss  of  generality,  that  N  >_  m  .  Theorems  4.1,  4.2  and  4.3 

characterize  the  limiting  d.f.'s  of  these  order  statistics  and  extend  the 
results  of  Berman  [2]  who  obtained  the  limiting  d.f.'s  of  the  maximum  under 
such  conditions.  Let: 

(4.1)  pjn_lNn  1  yj  -  An(y) 

and  assume  that  there  exists  a  d.f.  A(y)  ,  which  may  be  improper,  such  that 
for  all  y  in  its  continuity  set: 

(4.2)  litn  A^(y)  ■=  A(y)  . 

n-+°° 

To  ensure  that  the  limiting  d.f.'s  are  proper,  A(y)  will  be  required  to 
satisfy  one  of: 

(4.3)  A(0+)  -  A(O-)  -  0  ;  A(«)  -  A(O-)  >  0 

(4.4)  A(0+)  -  A(O-)  <  1  ;  A(»)  -  A(O-)  -=  1 

(4.5)  A (Of)  -  A(O-)  -  0  ;  A(»)  -  A(O-)  =  1  . 

Theorem  4.1: 

If  F(x)  is  in  the  domain  of  attraction  of  $,,..(x)  ,  then  -  b 

(1)  n  \  n  i 

has  a  limiting  d.f.  G(x)  which  is  necessarily  of  the  type: 


« 


t 


G(x)  ■=  0 


for  x  <  0 


00 

m-1  f .  a.j 

-  I  /V 

J-O-'n  J 


exp  [-yx  ]dA(y) 


for  x  >  0  ,  a  >  0 


where  A(y)  satisfies  condition  (4.3), 


Proof : 


Let  G  (a  x  +  b  )  =  P  la  ^  (^’n  -  b  )  <  x  ,  so  that 

nm  n  n  (  n  \  n  n  /  —  ) 


(4.6)  G 


(a  x  +  b  )  =  E  <  y  ^  (a  x  +  b  )F^  (a  x  +  b  ) 

i  n  n  n  n  n  n 


where  E^O)  denotes  expectation  with  respect  to  N  and  the  constants 

and  {b  }  are  such  that: 
n 


(4.7) 


lim  F(a  x  +  b  )  «  ^...(x)  ,  V  x  . 
n  n  (.17 


It  is  easy  to  show  that  lim  G  (a  x  +  b  )  *  0  for  x  <  0  :  consider  one  of 
1  nm  n  n  — 

the  terms  in  the  sum  (4.6)  with  x  >  0  : 

<4'8)  E„!("rj(V  +  bn)  +  bn>!  ■ 


1  f"F(V  +  V|3„  (n(N-  U...(N- 


Jl  (haDx  +  V) 


FN(a  x  +  b  ) 
n  n 


From  (4.7),  lim  F(a  x  +  b  )  =  1  and  lim  nF(a  x  +  b  )  »  -log1!1  N  (x)  ,  so 

n  n  n  n  6  (1) 

for  the  first  factor  in  (4.8) 


(4.9) 


nF(a  x  +  b  ) 
n  n 


n  n 


mry(  4  '•  l08'(i)(x))J 
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Now 


lim  ENjn"JN(N  -  1) .  . .  (N  -  j  +  l)e“Sn  Nj  *=  lim  |(n_1N)^  e_Sn  N| 


(4.10) 


oo 

lim  J'  y^e  SydAn(y) 


OO 

J'  y^e  SydA(y) 


for  s  >  0  ,  by  the  extended  llelly-Bray  lemma.  What  has  been  shown  so  far  is 


(4.11) 


00 


g(y)dAn(y) 


g(y)dA(y) 


where  g(y)  =  y^  exp[-sy]  ;  now  let  g^(y)  c  y^  exp[ny  logFfa^x  +  b^) ]  . 

Then  since  -n  log  ^(a^x  +  b^)  -*  -log$^(x)  >  0  ,  for  x  >  0  ,  it  may  be 
shown  that  g^(y)  **  8 (y )  uniformly  in  y  for  each  fixed  x  >  0  as  n  -*•  00  , 
with  s  =  -log<t'^j(x)  .  Hence 


(4.12) 


f 8n(y)dAn(y) 

*/o 


g (y)dA(y) 


which  is  the  limit  of  the  second  factor  in  (4.8);  combining  this  with  (4.6), 
(4.9),  (4.12)  and  (1.1),  the  main  part  of  the  theorem  follows. 

Uniqueness  of  the  limit  d.f.  follows  from  Khintchine's  theorem  [7]  that 
a  sequence  of  d.f.'s  {F^Ca^x  +  ^n) )  can  have  at  most  one  proper  limit  type 
(4.3)  follows  from  the  conditions  G(0+)  <  1  and  G(°°)  =  1  that  the  limit 


d.f.  be  proper. 
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Theorem  4.2: 

If  F(x)  is  in  the  domain  of  attraction  of  *  (x)  ,  then 

-1  /  (m)  .  \  (2)  ‘  ' 

°n  \  n  ”  bn/  bas  a  d.f.  G(x)  which  is  necessarily  of  the  type 

oo 

G(x)  "  1  "  1/ }  cxp[-y(-x)“a)dA(y)  for  x  <  0 

J“u  Q 

m  1 

for  x  ^  0  ,  a  >  0 

where  A(y)  satisfies  condition  (4.4). 

Theorem  4.3: 

If  F(x)  is  in  the  domain  of  attraction  of  <t>  (x)  ,  then 

-1/  (m)  \ 

fln  l^n  bn)  has  a  limitln8  d.f.  G(x)  which  is  necessarily  of  the  type 
G(x)  “  1  -  l  exp[-yeX]dA(y)  —  <  x  <  » 

where  A(y)  satisfies  condition  (4.5). 

The  proofs  of  Theorems  4.2  and  4.3  are  analogoos  to  that  of  Theorem  4.1. 
Conditions  (4.4)  and  (4.5)  follow  by  considering  C(O-),  G(-«)  and  G(— ), 

G(+«)  respectively. 
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5.  APPLICATION  TO  SUPERPOSITION  OF  RENEWAL  PROCESSES 

Necessary  and  sufficient  conditions  for  a  sum  of  a  large  number  of  re¬ 
newal  processes  to  be  close  to  a  (nonhomogeneous)  Poisson  process  have  been 
given  by  Grigelionis  [8],  he  considers  n  general  renewal  processes,  where 
Fnr(x)  and  Fnf(x)  are  respectively  the  d.f.'s  of  time  to  the  first  event 
and  between  subsequent  events,  1  <_  r  n  .  Under  the  condition 

lim  max  max  {F  (x)  ,  F  (x) }  ■=  0  ,  for  all  x  ,  he  proves 
n  r 

Theorem  5.1: 

For  the  superposed  process  to  converge  as  n  to  a  Poisson  process 

with  parameter  A(x)  ,  it  is  necessary  and  sufficient  that 

n  „ 

l  Fnr  ^x)  A(x)  V  x  >  0  . 

r-1 

It  is  of  interest  to  know  what  form  the  function  A(x)  may  take  when 
the  underlying  renewal  processes  are  identical  and 

F  (x)  ■  F  (x)  *  F(a  x+b)  ,  l<r<n,  where  {a  >  0}  and  {b  }  are 
nr  nr  n  n  —  —  n  n 

suitably  chosen  constants.  F(x)  is  a  fixed  d.f.  on  the  nonnegative  real 
axis,  and  it  is  convenient  to  assume  that  F(e)  >  0  for  e  >  0  and  take 
b  ■  0  .  The  following  result  is  suggested  by  Theorem  2.1. 

Theorem  5.2: 

If  the  sum  of  n  identical  renewal  processes  converges  to  a  Poisson 
process  with  parameter  A(x)  ,  then  A(x)  -  axa  for  some  a  >  0  ,  a  >  0  . 
Moreover,  F(x)  is  in  the  domain  of  attraction  of  $^(x)  • 
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Proof : 

Observe  that  the  first  event  in  the  sum  process  is  the  same  as  the 
minimum  of  a  set  of  n  independent  and  identically  distributed  random 
variables,  so  that  Gnedenko's  results  show  that  this  minimum  has  a  limiting 
d.f.  which  is  one  of  those  in  (1.1).  However,  the  d.f.  ®(2)^X^  Can 
ruled  out  because  its  domain  of  attraction  contains  only  d.f.'s  corresponding 
to  random  variables  unbounded  below  (6).  Moreover,  8  d.f.  on 

the  whole  real  axis  and  so  cannot  generate  a  renewal  process.  Thus  the 
limiting  d.f.  of  the  minimum  must  have  a  d.f.  of  the  same  type  as 
which  corresponds  to  a  Poisson  process  with  parameter  A(x)  =  axa  . 

In  fact,  the  same  form  for  A(x)  holds  under  wider  conditions  than 
those  above.  Suppose  F(x)  is  in  the  domain  of  attraction  of  4>^(x)  — see 
Section  1 —  so  that 


(5.1) 


lim  F(tx)  /  F(x)  =  t  ,  a  >  0 
x>0+ 


Let 

(5.2)  F  (x)  =  F  (x)  «=  F(y  a  x)  , 

nr  nr  '  r  n  * 

so  that  in  the  terminology  of  Section  2  one  may  think  of  identical  components 
which  are  subjected  to  different  levels  of  stress  in  different  sockets. 
Assuming  that 


(5.3) 

then  if 


£  ya  "  00  and  y 
r-1  r  n 


aR  are  chosen  to  satisfy 
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(3.4) 


n  ' 

an  =  ^inf  x  |  F(x)  >_  lj\  ^  y“ 


It  may  be  shown  (see  [10])  that 


(3.5; 


7  F(y  a  x)  -*■  x 
L,  r  v> 
r=l 


Hence,  Grlgelionis'  conditions  are  satisfied,  and  one  has: 

Theorem  5.3: 

The  sum  of  n  renewal  process  of  the  form  specified  by  (5.1),  (5.2) 
and  (5.3)  converges  to  a  Poisson  process  with  parameter  A(x)  =  axa  . 

Now  suppose  that  the  number  of  renewal  processes  is  itself  a  random  vari¬ 
able  N  •  *■  the  case  of  order  statistics  in  Section  4.  Using  the 

n 

notation  of  that  section: 


Theorem  5.4: 

The  sum  of  a  random  number  of  identical  renewal  process  of  the  form 

F  (x)  ■  F  (x)  *=  F(a  x)  ,  where  F(x)  satisfies  (5.1)  and  {a  >  0}  are 
nr  nr  n  n 

suitable  constants,  converges  to  a  counting  process  as  n  -*■  00  in  which  the 
d.f.  Gm(x)  of  the  time  to  the  m1^  event  is  given  by 

G  (x)  0  for  x  <  0 

m  — 

m-l  f ,  ouj 

■  1  -  \  I  '^7,  exp[-  yxa]dA(y)  for  x  >  0 

J-0,b 

where  A(y)  ■  lim  pjn  <_  yj  satisfies 

(5.5)  A  (Of)  -  A(O-)  -=  0  ;  A(«)  -  A(O-)  ■=  1  . 


r 
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Proof : 

A  corresponding  result  for  order  statistics  from  a  random  sample  has 
been  proved  in  Theorem  4.1;  it  remains  to  be  shown  that  the  d.f.  of  the  m**1 
order  statistic  in  such  a  random  sample  is  asymptotically  the  some  as  the 
d.f.  of  the  mt^1  event  in  the  superposed  process.  From  (5.5),  k  and  K  may 
be  chosen,  0<k<K<°°,so  that  pjk  <_  n  <_  k|  >  1  -  c  ,  where  c  >  0 
is  arbitrary.  Then  for  each  in  this  intt  rval  the  asymptotic  closeness 

of  the  two  limiting  d.f.’s  in  question  may  be  shown  as  in  the  proof  of 
Theorem  2.3.  Unconditioning  on  and  letting  c  0  ,  the  desired  result 

is  obtained. 

Theorem  5.4  shows  that  randomness  of  the  number  of  renewal  processes 
leads  to  counting  processes  which  are  mixtures  of  (nonhomogeneous)  Poisson 
processes. 

Example: 

Take  a  ■  2  ,  A(y)  -  1  -  e  ^  .  The  limiting  d.f.'s  are 

G  (x)  ■=  0  for  x  <  0 

m  “ 

■  {x2/(l  +  x2)}m  for  x  >  0 


33 


6.  COMPOSITION  Ob'  COHERENT  STRUCTURES 

Many  multi-component  systems  occuring  in  the  field  of  reliability  are 
coherent  structures  (e.g. ,  see  [1]).  A  coherent  structure  consisting  of  .n 
components  ci,c2’  •  •  •  »  cn  is  described  by  its  structure  function  <j>(x) 
which  is  1  if  the  system  works  and  0  otherwise,  where  x  *  ....  x^) 

and  is  1  if  component  c^  is  working  and  0  if  it  is  failed.  The  binary 

function  <p  is  required  to  satisfy 

(i)  Each  component  c^  must  be  e88ential\  that  is  to  say,  there  is 


a  realization  x^,  .. 

•’  xi-l  ’ 

*1+1’  •" 

such  that 

4>(xlf  •  •  •  ,  ,  1  , 

xi+l*  ••• 

,  xn)  -  1 

and 

$ (Xj,  • • • ,  ®  » 

xi+r  ••• 

•  V  '  0 

• 

(ii)  <f>(x)  1  4>(y)  for  all 

X  =  (Xj, 

....  xn) 

and 

y  “  . yn> 

such  that  x^  _<  y^  ,  1  <_  i  <_  n  . 

From  (i)  and  (ii)  it  follows  immediately  that  $(1)  «=  1  and  ip(0)  =0  , 
where  1  -  (1,1,  ....  1)  and  0  =  (0,0,  ...,  0)  . 

Most  of  the  classical  problems  in  extreme  value  theory  are  contained  in 
the  following  general  question:  "Given  a  coherent  structure  with  a  finite 
number  of  components  and  some  procedure  to  increase  the  number  of  components 
without  bound,  what  are  the  possible  limiting  d.f.'s  for  the  system  lifetime 
for  given  component  lifetime  d.f.'s?"  For  example,  the  minimum  of  a  set  of 
random  variables  corresponds  to  the  lifetime  of  a  series  structure,  and  the 
maximum  likewise  corresponds  to  a  parallel  structure.  The  so-called  k-out-of-n 
structure  generate  order  statistics  an'1  series-parallel  and  parallel-series 
systems  give  rise  to  the  minimax  and  raximin  (see  [3])  of  double  arrays  of 
random  variables. 

In  this  section,  results  are  obtained  when  the  method  of  "expanding"  the 
structure  is  that  of  repeated  composition.  Thus,  starting  with  a  system  of 
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n  components,  each  component  is  replaced  by  a  similar  system  giving  a  sy^  tern 

2  I 

of  n  components,  then  each  component  in  this  system  is  replaced  by  a  replica 

3 

of  the  original  system,  giving  a  system  of  n  components,  etc.  The  first 
three  stages  of  composition  for  a  three  component  system  are  illustrated  in 
Figure  6.1. 

It  is  assumed  that  all  components  are  independent  and  identical.  In 
particular,  let  p  *=  P{X^  *=  1}  be  the  reliability  of  each  component  c^  , 
where  is  the  binary  random  variable  designating  the  state  of  component 

c^  .  The  structure  function  $(*)  becomes  a  binary  random  variable  and  the 
reliability  of  the  structure  will  be  denoted  by  h(p)  *■  P { (X)  =  1}  .  Now 
suppose  that  the  components  fail  over  time  and  have  a  common  lifetime  d.f. 

F(t)  with  binary  indicator  random  variables  X^(t)  which  are  1  or  0 
according  as  component  c^  is  in  a  working  or  failed  state  respectively  at 
time  t  .  Then  P{X^(t)  *=  1}  *=  F(t)  and  so  the  system  has  lifetime  d.f. 

1  -  h(F(t))  .  The  system  which  results  after  n  compositions  has  lifetime 
d.f.  1  -  h^(F(t))  ,  where  h^  denotes  the  n-fold  composition  of  the 
function  h  with  itself.  It  is  assumed  throughout  this  section  that 
h(p)  t  p  .  Suppose  that  for  suitable  {a^  >  0}  and  {b^} 

(6.1)  h(n)(F(a  x  +  b  ))  -  G (x)  . 

n  n 

Then 

Lemma  6.1: 

If  G(x)  is  a  limiting  d.f.,  then  h(G(x))  =  G(ax  +  0)  for  some 
a  >  0  ,  6  . 

Remark : 

This  lemma  generalizes  the  concept  of  maximum  stable  d.f.'s  introduced 

by  Gnedenko.  The  relation  h^(G(x))  *  G(a  x  +  6  )  for  some  a  >  0  .  6 

n  n  n  ’  n 


I 


n  ■  3 


Figure  6.1 
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nnd  for  all  n  is  no  more  general,  since  it  is  implied  by  the  lemma  with 

a  ■  an  and  0  =  (1  +  a  +  . . •  +  aU  ^)6  . 

n  n 

Proof: 

From  (6.1), 

h(n+1)(F(a  x  +  b  ))  ->  h(G(x)) 
n  n 

and 

h(n+1)(F(an+1x  +  bn+1))  -*  G (x)  . 


But  by  Khintchine's  theorem  [7],  these  two  d.f.'s  must  be  of  the  same  type, 
i.e.  , 


h(G(x))  *  G(ax  +  B)  . 


If  a  ^  1  ,  set  G  (x)  «■  G 


f +  ■ 


so  that 


(6.2)  h(G  (x))  -  h^G^x  +  ^  "a))  “  G  |ax  +  ^  ^  ^  j  -  G  (ax)  . 

ft  ft  o 

If  a  -  1  ,  set  G  (x)  ■  G(logx)  and  a  *»  e  ,  so  that 

(6.3)  h(G  (x))  -  h(G(log  x))  »  G(log  x  +  log  a  )  -  G  (a  x)  . 


Theorem  6.1: 

G(x)  is  a  limiting  d.f.  if  and  only  if  G(x)  is  type  equivalent  to 
<Kx)  or  <Ji(e  )  ,  where  1  -  \p(x)  is  a  d.f.  and 


(6. A) 


h(^(x) )  ■  '^(ax) 


for  all  x  ,  where  a  >  0  ,  a  j*  1  . 
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Proof : 

t  /a  j'i  (6  3)  and  Lemma  6.1.  F°r  sufficiency, 
The  necessity  part  follows  from  (6.2),  (6.3) 

.  .  .  f  is  m  its  own  domain  of  attraction  by  choosing, 

notice  that  each  such  d.f.  is  m  its 


a""  ,  b  -  0  or  an 


1  b  -  -n  log  a  respectively. 

*  n 


The  following  1»  *  well-known  result  (e.C.  see  [U,  P-  199) 


Theorem  6.2: 

Excluding  the 
one  of  three  kinds 


trivia!  case  h(p)  -  p  .  the  reliability  function  h(p)  is 


(i) 

h(p) 

<  P 

for 

0  <  p  <  1 

(ii) 

h(p) 

>  P 

for 

0  <  p  <  1 

(iii) 

h(p) 

<  P 

for 

0  <  p  <  P0 

“  P 

Ko 

for 

P  ■  P0  is 

>  P 

for 

po  r  p  r  1 

Defining  h*(p>  -  1  -  «1  -  ,)  -  ♦*<*>  -  1  - 

h(p)  <  p  <*>  h*(p)  >  P  «nd 

(6.5)  hV<*»  -  1  -  l'«(-x))  -  1  -  ♦  <-»>  ‘  ’ 

so  that  it  suffices  to  consider  cases  (1)  and  (ill)  only. 


Lemma  6.2: 

If  h(p)  <  P  *nd 
situations  can  occur 


)  satisfies  (6.4),  then  neither  of  the  following 


(i)  a  >  1  ,  x  <  0  ,  0  <  *00  <  1 

(ii)  a  <  1  ,  x  l  0  ,  0  <  *00  <  1 
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Proof : 

(i)  i>(ax)  >_  ^(x)  ,  since  ax  x  and  i|>  is  monotone,  and 

hi (4/ (x ) )  <  \p(x)  by  the  assumptions.  Thus  with  (6. A),  a  contradic¬ 
tion  results.  The  proof  of  (ii)  is  similar.  Likewise,  one  can 
show 

Lemma  6.3: 

If  h(p)  is  S-shaped  about  p^  and  i»(x)  satisfies  (6.4),  then  none 
of  the  following  can  hold 

(i)  a  <  1  ,  x  >_  0  ,  0  <  ip(x)  <  p^  . 

(ii)  a  <  1  ,  x  <  0  ,  p  <  i*/(x)  <  1  . 

—  o 

(iii)  a  >  1  ,  x  _<  0  ,  0  <  ip  (x)  <  p^  . 

(iv)  a>l,x_>0,p  <  >i(x)  <  1  . 

From  (i),  (ii)  and  the  fact  that  ijj(x)  is  the  tail  of  a  proper  d.f., 
it  is  seen  that  (6.4)  cannot  hold  with  a  <  1  when  h(p)  is  S-shaped.  Thus 
the  solutions  to  (6.4)  take  one  of  the  following  five  forms 

la)  h(p)  <  p  ,  a  >  1  ,  \p  (x)  «=  1  for  x  <  0 

0  <  \p  (x)  <  1  for  x  >  0  . 

lb)  h(p)  <  p  ,  a  <  1  ,  0  <  ^ (x)  <  1  for  x  <  0 

tp(x)  ■  0  for  x  ^  0  . 

2a)  h(p)  S-shaped  about  p  ,  a  >  1  ,  'i'(x)  “  1  for  x  0 

0  <  <J<(x)  <  PQ  for  x  >  0 

2b)  h(p)  S-shaped  about  pQ  ,  a  >  1  ,  p&  <  ij/(x)  <  1  for  x  <  0 

<Kx)  -  0  for  x>0 

2c)  h(p)  S-shaped  about  pQ  ,  a  >  1  ,  pQ  <  ij/(x)  <  1  for  x  <  0 

0  <  ii(x)  <  P0 


for  x  >  0 
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The  class  of  all  solutions  to  (6.4)  can  now  be  generated;  for  example: 


Theorem  6.3: 

The  class  of  all  solutions  of  the  form  la)  may  be  obtained  by  the 
following  procedure: 

(i)  Pick  any  a  >  1  . 

(ii)  Pick  any  x^  >  G  . 

(iii)  Define  U/(x)  in  the  half-open,  half-closed  interval  (x  ,ax  ] 

o  o 

so  that  it  is  nonincreasing  and  not  constant  over  the  interval. 

Extend  ip (x)  to  the  positive  real  axis  by  defining 

(k)  (  /  -k  \)  k  k+1 

<|)(x)  ■  h  j^a  x  jj  ,  for  a  xq  <  x  <_  a  XQ  »  where  a  negative 

index  k  means  the  (-k)*'*1  composition  of  the  inverse  function 


Similar  procedures  may  be  derived  for  solutions  which  are  of  the  other 
forms . 

It  is  clear  that  the  class  of  limit  d.f.'s  thus  generated  is  very  large. 
In  the  case  of  composition  of  series  structures,  the  class  of  limit  d.f.'s 
has  been  characterized  by  Mejzler  [13]  and  includes  d.f.'s  other  than  those 
of  the  type  (1.1).  For  example,  the  function 


(6.6) 

X 

H 

O 

for  x  <_  0 

3 

■  exp[-  x  (5  +  sin(2ir  log  x'log  2))] 

for  x  >  0 

g 

is  the  tail  of  a  d.f.  and  satisfies  (6.4)  for  h(p)  *=  p  ,  a  *  2  . 
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